Composition of fluctuations of different observables 
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We demonstrate that description of fluctuations observed in multiparti- 
cle production processes using Tsallis statistics approach (in which fluctu- 
ations are described by the nonextcnsivity parameter q) leads to a specific 
sum rule for parameters q seen in different observables which can be verified 
experimentally. 
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1. Introduction 

When analyzing multiparticle production data the standard tool used is 
statistical modelling [1]. However, this approach does not account for the 
possible intrinsic nonstatistical fluctuations in the hadronizing system which 
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usually result in a characteristic power-like behavior of the single particle 
spectra and in the broadening of the corresponding multiplicity distributions 
(and which can signal a possible phase transition(s) [2]). One possibility 
to include this and still remain in the domain of a statistical approach 
is to use the so called Tsallis statistics O HI [5] (represented by Tsallis 
distribution, hq{E)) which accounts for such situations by introducing in 
addition to the temperature T, a new parameter, q> 1, directly connected 
to fluctuations (SJ [7] (for g — )• 1 one recovers the usual Boltzmann-Gibbs 
distribution, f{E)): 
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^ie.p(-f), (2, 

The most recent applications of this approach come from the PHENIX Col- 
laboration at RHIC [8] and from the CMS Collaboration at LHC [9]. The 
parameter q is entirely given by intrinsic fluctuations in the system, cf. Eq. 
(gl) helovB- 

Before proceeding any further we must emphasize two points. First, the 
relation between parameter q and fluctuation of temperature was derived 
in [6], where it was shown that starting from some simple diffusion picture 
of temperature equalization in a nonhomogeneous heat bath (in which local 
T fluctuates from point to point around some equilibrium temperature, T) 
one obtains an evolution of T in the form Langevin stochastic equation and 
distribution of 1/T, /(1/T), as solution of the corresponding Fokker-Planck 
equation. It turns out that /(1/T) has the form of a gamma distribution, 



1 \q-lTj \ q-lT 



where o — 1 = — -. 

(i/r)2 



^ See also [10] . One must admit at this point that this approach is subjected to a 
rather hot debate of whether it is consistent with equihbrium thermodynamics or 
else it is only a handy way to phenomenologically describe some intrinsic fluctuations 
in the system |Tl]. However, as was recently demonstrated on general grounds in [12| . 
fluctuation phenomena can be incorporated into a traditional presentation of thermo- 
dynamic and Tsallis distribution, Eq. belongs to the class of general admissible 
distributions which satisfy thermodynamical consistency conditions and which are 
therefore a natural extension of the usual Boltzman-Gibbs canonical distribution Eq. 
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Convoluting exp(— i?/T) with such a /(1/T) one obtains immediately Tsal- 
hs distribution hq{E) ^ The parameter ^, i.e., according to Eq. ([H 
also the temperature fluctuation pattern, is therefore fully given by the pa- 
rameters describing this basic diffusion process (cf., [6] for details)!. 

The second point is that, as was shown in [15], temperature fluctuations 
in the form given by Eq. ([3]) result in an automatic broadening of the 
corresponding multiplicity distributions, P{N), from the poissonian form 
for the exponential distributions, Eq. ([2]), 

(N)^ - - E 

P(N) = exp {-N) where N = —. (5) 

to the negative binomial (NB) form for the Tsallis distributions, Eq. ([1]) 
(cf., [15], for details), 

P(Ar) = + ^ ' where k = (6) 

^ ^ r{N+l)T{k) l^^^{Niy^+'')' q-l 

Notice that in the limiting cases of g — >• 1 one has A; — t- oo and ([6]) becomes 
a poissonian distribution whereas for g — >■ 2 on has /c — )• 1 and ([6]) 
becomes a geometrical distribution. It is easy to show that for large values 
of N and (N) one obtains from Eq. ^ its scaling form, 

/ N \ k^ 
{N)P{N) ^i;^z = ^j= ^/-i eM-kz), (7) 

in which one recognizes a particular expression of Koba-Nielsen-Olesen (KNO) 
scaling [16j 0. 



2. Results 

Proceed to a detailed analysis of different observables, like P{N), dN/dy 
or dN/dpT, fluctuations in which are expected to differ from each other 
and therefore to result in different values of the corresponding parameters 

This was recently generalized to account for the possibility of transferring energy 
from/to a heat bath, which appears to be important for AA applications [HITSj and 
for cosmic ray physics however, we shall not discuss this issue here). 
^ It is worth mentioning at this point that, as shown in [T7], fluctuations of N in the 
poissonian distribution (O taken in the form of tp{N / < N >), Eq. (Q, lead to the 
NB distribution ©. 



4 



q'sum'rule printed on February 23, 2011 



q. Indeed, from our experience with pp collisions [18j we know that one 
can obtain a very good description of the whole range of pT spectra (oc 
expg {-pt/T) with (T = Tt [GeV]; qr) = (0.134; 1.095), (0.135; 1.105) and 
(0.14; 1.11) for energies (in GeV) 200, 540 and 900, respectively 0. These 
values should be compared with the corresponding values oi {T = T^^q = 
qi) obtained when fitting rapidity distributions (oc exp^ (— /i^ cosh y/T)) 
at the same energies: (11.74; 1.2), (20.39; 1.26) and (30.79; 1.29). It was 
noticed there that — 1 has the same energy behavior as \/k in the NB 
distribution fitting the multiplicity distributions at corresponding energies 
{qL — 1 = —0.104 + 0.058 In \/i). This means that fluctuations of the total 
energy are in this case mainly driven by fluctuations in longitudinal phase 
space. An explanation proposed in ^18j was the following. Noticing that 
q — 1 = a'^{T)/T'^ (i.e., it is given by fluctuations of total temperature T) 
and assuming that cr'^iT) = a'^{TL) + a'^{TT), one can estimate that the 
resulting values of q should not be too different from 



The situation is noticeably different for nuclear collisions, which we shall 
now addres^. As shown in Fig. [H data for dN/ dy and dN/ dfix can be fitted 
perfectly by means of Eq. ([1]). However, the behavior of the q parameters 
obtained is quite interesting, as displayed in Fig. [2J At first, the parameter 
q from P{N) turns out to depend on the centrality of the collision defined 
by the number of participants of projectile, Np (left panel of Fig. [2]), 



{A - mass number of colliding nuclei and a = 0.98) [13J. Whereas for small 
centralities it approaches situation encountered in pp, the more central the 
event, the smaller is (7 — 1, i.e., the nearer to a poissonian the corresponding 
P(N). This time for each centrality qx are larger than q (both for results 
based on NA49 data [20] and from fits presented in [22] based on PHENIX 
data |21]). In right panel of Fig. [2] we collected all results for the most cen- 
tral events from NA49 [20j. Notice that they clearly display opposite trend 

* There is recent compilation of essentially all results for pr spectra, including recent 
LHC data [19]. It shows that q{s) — 1.25 — O-SSs"""^"*, which nicely reproduces 
results mentioned here. 

^ For this purpose we use mainly NA49 data on Pb + Ph collisions [20] because, at the 
moment, only this experiment measures both (at least for the most central collisions) 
multiplicity distributions, P{N), and distributions in rapidity y, transverse momenta, 

Pt, and transverse masses, /it = (jn^ + p^^ and this property is crucial for further 
considerations. PHENIX results [3T] analyzed in [52] are also shown for comparison. 
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Fig. 1. (Color online) Fits using Eq. ((T)) to ^20] data for dN/d^T (a) and for dN/dy 
(b) for central collision Pb + Pb at different energies. 



to that encountered for the pp collisions mentioned above: both, and qt 
(obtained from pT distributions are now greater than q and have (approxi- 
mately) a visible similar dependence on Np. However now q^ < qj-, again, 
this is opposite to what was seen in pp. 

A natural question is, what causes such different behavior of the parame- 
ter q in this case? The answer we propose: When extracting values of q from 
the rapidity distributions a tacit assumption was that fix in E = /xy cosh y 
remains constant (i.e., it does not fluctuate). However, this is too crude, 
because data show that fluctuates as well. To account for this fact notice 
that expq (—E/T) = exp^ [— {^t/T) coshy] = expg(— zcoshy), i.e., that fits 
to rapidity distributions provide us with fluctuations not so much of parti- 
tion temperature T but rather of the variable z = ^ix/T. This in turn can 
be written approximately as: 

Var{z) -jl^^Var M + ^ • (10) 
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Fig. 2. (Color online) (a): q for different centralities measured by the number of 
projectile participants, Np. Here q from P{N) are our results from Var{N)/ (N) 
[l3] to be compared with q — qx obtained from dN/ dpr data of NA49 (cf. first work 
of [20]). PHENIX results [Slj analyzed in [22] are also shown for comparison. (6): 
Energy dependence of q for the most central events. All results were obtained for 
the sake of this presentation using distributions provided by [2D], i.e., respectively, 
dN/dfiT, dN/dy and dN/dpT- The errors are similar to those presented as an 
example for q = Ql obtained from dN/dy. Open symbols correspond to uncorrected 
values of q, full symbols to values corrected by means of the procedure proposed 
in this work and explained in the text. 



Because 



and because 



\ / (^j^^ (1/T)2 (T)2 ^ ^ 



Var{z) _ Varifix) , Var(T) 

+ /rr\9. U^j 
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one obtains following sum rule connecting different fluctuations 

_ dg Var{T) ^ Var{z) _ Var (fir) , . 

This sum rule is our main result and its action is presented in the right 
panel of Fig. [21 It connects total q, which can be obtained from the 
analysis of the NB form of the measured multiplicity distributions, P(N), 
with qL — 1 = Var{z)/ (z)'^ , obtained from fitting rapidity distributions 
and Var^HT) / (iJ't)'^ obtained from data on transverse mass distributions. 
When extracting q from distributions of dN/d^T we proceed analogously 
with z = coshy/r. 





o p+p at s'°= 540 GeV 




Tln= 35.1, q= 1.483 
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Fig. 3. (Color online) (a): Example of best fit to dN/dy for Ecu = 540 GeV 
(data are the same as in [IZ]). {h): The q-dependence of the admissible T/ ^t- (c): 
^'-dependence of the corresponding . The full red dot shows the minimal value 
of x^i see text for details. 



Some explanatory remarks on the apparent discrepancies between pp 
and AA data are in order here. When fitting dN/dy data on pp in [T7] ht 
was kept constant (and given by the {px) for given energy). This means 
that all effects related to its fluctuations was attributed to the fluctuations 
of the "partition temperature" T = Tl (which is therefore only one of the 
fitted parameters and, for example, it cannot be used to calculate the mean 
energy). This was possible because, as observed in [SS^jfS], in the fitting pro- 
cedure parameters T and q are strongly correlated. To illustrate this fact 
we present in Fig. Elan example of the best fit to dN/dy for Ecm = 540 
GeV together with g-dependency of T / iit and parameter representing 
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the goodness of the fil|f|. Notice that for different values of q and T/ht 
we can describe dN/dy with reasonable accuracy. Parameters T and q are 
strongly correlated, for Eqm = 540 GeV shown here T{q)/ fix = 56.5—14.5(7. 
Comparing values of T/jit and q shown in Fig. [3] to those reported for this 
energy earlier [T7] {q' = 1.26 and T'/fiT = 35.1) one can observe that the 
difference Aq = q — q' = 0.22 is roughly the same as the correction caused 
by fluctuations of discussed abov^. 



3. Summary 

We have demonstrated that fluctuations of temperature T, together with 
fluctuations of other variables, result in the sum rule formula, Eq. ([T3|), con- 
necting q obtained from an analysis of different distributions. This allows 
us to understand why in AA collisions fluctuations observed in multiplicity 
distributions are much smaller than the corresponding ones seen in the ra- 
pidity distribution or in the distribution of transverse momenta (i.e., why 
the corresponding q parameters evaluated from distributions of different ob- 
servables are different). This issue should be checked further when complete 
sets of data become available from the experiments at LHC (especially from 
ALICE). 
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